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i-^ ' Abstract. In the CDM scenario, dark matter halos are assembled hierarchically from smaller 

subunits. A long-standing problem with this picture is that the number of sub-halos predicted 
^ , by CDM simulations is orders of magnitudes higher than the known number of satellite galaxies 

^ ' in the vicinity of the Milky Way. A plausible way out of this problem could be that the majority 

I of these sub-halos somehow have so far evaded detection. If such "dark galaxies" do indeed exist, 

. gravitational lensing may offer one of the most promising ways to detect them. Dark matter 

sub-halos in the 10*^ - 10^° Mq mass range should cause strong gravitational lensing on (sub- 
)milliarcsecond scales. We study the feasibility of a strong lensing detection of dark sub-halos 
by deriving the image separations expected for density profiles favoured by recent simulations 
and comparing these to the angular resolution of both existing and upcoming observational 
facilities. We find that there is a reasonable probability to detect sub-halo lensing efi^ects in 
^ I ^ ' high resolution observations at radio wavelengths, such as produced by the upcoming VSOP-2 

\l I satellite, and thereby test the existence of dark galaxies. 
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1. Introduction 

A multitude of recent cosmological studies suggest that most of the matter in the universe 
consists of non-baryonic particles which are believed to constitute the cold dark matter (CDM) 
' [T]. While the CDM scenario has been very successful in explaining the formation of large-scale 

structures in the Universe (see e.g., [2], for a review), its predictions on the scales of individual 
galaxies have not yet been confirmed in any convincing way. One particularly interesting feature 
of the CDM model is the high level of halo substructure generated, seemingly inconsistent with 
the known number of satellite galaxies in the vicinity of the Milky Way (e.g., [3', '5]). This 
discrepancy persists even when considering the newly discovered ultra-faint dwarf galaxies in 
the Sloan Digital Sky Survey and correcting for its sky coverage [6] . One possible way out of this 
problem is to assume that most of these low-mass halos correspond to so-called dark galaxies [7], 
i.e., objects of dwarf-galaxy mass which either do not contain baryons or in which the baryons 
have not formed many stars. 

Gravitational lensing may in principle offer a route to detecting even completely dark galaxies. 
Dark matter sub-halos in the 10^-10^^ Mq range are expected to cause gravitational milli-lensing, 
i.e., image splitting at a characteristic separation of milliarcseconds [8l[9]. To put the CDM sub- 
halo predictions to the test, it has been suggested that one should target quasars which are 



already known to be gravitationally lensed on arc second scales, as one can then be sure that 
there is a massive halo well-aligned with the line of sight, which substantially increases the 
probability for sub-halo milli-lensing [9] . Indeed, the magnification associated with milli-lensing 
has long been suspected to be the cause of the flux ratio anomalies seen in such systems [10, llj. 
Sub-halo milli-lensing has also been advocated as an explanation for strange bending angles of 
radio jets [12] and image positions which smooth halo models seem unable to account for [13]. 

Here, we take a critical look at the prospects for strong-lensing detections of dark sub-halos 
in the dwarf-galaxy mass range. In ^ we derive the image separations expected for sub-halo 
density profiles favoured by recent simulations. In ^ we estimate the overall sub-halo lensing 
probabilities for point sources as well as extended sources. Finally, we discuss several effects 
that could affect our predictions and present our conclusions in ^ 

Throughout the paper, we assume a ACDM cosmology with f^A = 0.762, 0,m = 0.238 and 
h = 0.73 {Hq = lOO^kms^^ Mpc^^) in concordance with the WMAP 3-year data release |14] . 

2. Image separations 

To first order, the image separation produced through strong lensing by an extended object with 
a density that decreases as a function of distance from the centre is given by 



where D^i represents the angular-size distance between observer and lens, and i?E represents the 
linear Einstein radius. The latter is defined as the radius inside which the mean surface mass 
density S of the lens equals the critical surface mass density 



where -Dos and D\s are the angular-size distances between observer and source, and lens and 
source, respectively. 

The resulting image separations will thus heavily depend on the sub-halo surface mass density 
profiles S(r). Previously proposed strategies to detect image splitting [HI [15] assume that the 
image separations of sub-halo lenses are similar to those produced by a singular isothermal 
sphere (SIS). Unfortunately, this assumption is difficult to justify since theoretical arguments, 
simulations, and observations do not favour this form of density profile for dark matter halos in 
the relevant mass range. We have studied the feasibility of strong-lensing detection of dark sub- 
halos by deriving the image separations expected for (more realistic) density profiles favoured 
by recent simulations. These are NFW [IB], M99 [IT], N04 [IB], H03 [TO] and K04 [20j where we 
also have taken the impact of stripping and truncation into account. Details on these density 
profiles and a comparison of the resulting mean surface mass density profiles S(< r), can be 
found in [21J. 

In Fig. dl we plot the image separations predicted for 10^-10^^ Mq sub-halos at a redshift 
z\ = 0.5 and a source at = 2.0. We compare this to the angular resolution of a number of 
planned or existing observational facilities, operating at a wide range of wavelengths. These 
include the proposed MAXIM pathfindei0 in X-rays; VLTI with the proposed VSI instrument 
[22], the planned and SIM PlanetQuesT!^ satellites in the optical/near-infrared; the 

^ maxim.gsfc.nasa.gov/pathfinder.html 

^ www.rssd.esa.int/index.php?project=GAIA 

^ planetquest.jpLnasa.gov/SIM/ 







Figure 1. Sub-halo mass versus image separation for those density profiles that give rise to 
image separations on scales larger than micro-arc seconds. The angular resolution of a number 
of existing and planned observational facilities have been indicated by horizontal dashed lines, 
marked with labels (see main text for details), a) The different diagonal lines represent SIS 
(thin solid), truncated SIS (thick solid), M99 (medium dashed), stripped M99 (thin dashed) and 
truncated M99 (thick dashed) sub-halo models, b) The different diagonal lines represent NFW 
(thin dash-dotted), truncated NFW (thick dash-dotted) and truncated N04 (solid with circles) 
sub-halo models. 



currently available EVPvEl, HSA|, VLBAI arraysplus the planned ALMAIII, EVLAI, SKAI arrays, 
and also space- VLBI with the planned VSOP-^^ programme at radio wavelengths. Please note 
that here we consider only the best resolution limits attainable with these telescopes, whereas 
the resolution at the wavelengths that maximise the number of observable high-redshift sources 
may be considerably worse. 

From Fig. [T] it becomes obvious that there are large differences between the image separation 
predictions of the various halo models. As the discrepancy between the number densities of 
luminous galaxies and dark matter halos does not start to become severe until the halo mass 
drops below 10^° M© (e.g., [3 [23]), sub-halos at masses below this limit need to produce 
measurable image separations (^ ^ 4 x 10~^" for VSOP-2, which has the best theoretical 
resolution among the telescopes included in Fig. [T|) in order for dark galaxies to be detectable 
through image-splitting effects. Out of the halo models tested, only two actually meet this 
criterion without adhering to sharp truncations: the SIS and the M99 halos. The II03 and K04 
profiles both give image separations smaller than 10~^" for all the halo masses considered and 
are therefore completely outside the plotted region. Even in the optimistic case of an M99 halo 
(in either its original or stripped form, whereas the sharp truncation is not considered realistic, 
see |21j for details), the image separations are a factor of ~ 3-7 smaller than those predicted 
for a SIS (and ~ 30-60 times smaller than those of a truncated SIS), rendering only the few 
most massive sub-halos (~ 10^*^ Mq or slightly higher) detectable at ~ 0.01" resolution (GAIA, 

www.evlbi.org 
^ www.nrao.edu/HSA/ 
^ www.vlba.nrao.edu 
^ www.alma.info 
* www.aoc.nrao.edu/evla 
® www.skatelescope.org 

www.vsop.isas.ac.jp/vsop2 



SIM, and ALMA). At milliarcsecond resolution (VLTI and SKA), dark galaxies with masses 
> 10^ Mq may become detectable. To probe further down the sub-halo mass function, sub- 
milharcsecond-resolution facilities {MAXIM pathfinder, HSA, EVN, VLBA or VSOP-2) will be 
required. 

These estimates are based on the assumption that the sub-halos can be treated as isolated 
objects. The effects of external convergence k and shear 7 due to the galaxy (and halo) hosting 
these sub-halos have thereby been neglected. The question remains whether the boost factor 
/boost = A0K,7/A^K=o,7=o from the host potential can become considerably large for some of the 
more realistic density profiles considered here. To investigate this, we use ray-tracing simulations 
to numerically assess the distribution of /boost- We concentrate on two sub-halo models that 
span the range from effectively undetectable (H03) to favourable for detection (M99). In both 
cases, the tail of the /boost distribution extends up to very high values, which means that it is in 
principle possible to find macro- images for which sub-halos of masses much lower than indicated 
by Fig. [1] can be detected through image-splitting effects. However, such macro-images are 
exceedingly rare. For sub-halos with M^uh = lO^'^ M©, the probability of having /boost ^ 100 is 
~ 10~^ in the case of M99 and « 0.03 in the case of HQS. The expectation value for the boost 
factor is (/boost) ~ 2.3 for the M99 sub-halo and (/boost) ~ 14 for the H03 sub-halo. This is 
enough to shift the A9 prediction for the M99 sub-halos fairly close to that of an un-truncated 
SIS in Fig. [H but insufficient to move the II03 prediction (A^k=o,7=o ~ 10~^^" at this mass) 
into the detectable range. 



3. Lensing probabilities 

We assume a sub-halo population following the model proposed in ^24j. The sub-halo abundance 
per unit host halo mass (ignoring the high mass cut-off, Mgub > O.OlMhost) can be evaluated by 



-1.9 

dm yh^^MQ^ 

In [25], this finding was confirmed and extended to lower sub-halo masses using a high 
resolution simulation of CDM substructure in a Milky Way-sized halo. Within a host halo there 
is then a fraction /sub(-^host) -^sub_min) -^sub_max) of its mass Mhost in the form of sub-halos with 
minimum and maximum masses Msub_min and Msub_max! respectively. Here we assume sub-halo 
masses Mgub in the range 4 x 10^ < MgubZ-^^^o < lO^'', corresponding to the interval probed in 
[25]. The total mass fraction of a host halo with M = 1.8 x IO^^Mq in sub-halos is then about 
5% for this sub-halo mass range. 

Within this simulation, it has also been shown that the sub-halo number density profile can 
be fitted by the following form [26j : 

n(< x) 12x^ 

N " 1 + 11x2' ^ ' 

where x is the distance to the host center in units of -R20O1 the radius at which the mean 
enclosed density equals 200 times the mean mass density of the Universe, n(< x) is the number 
of sub-halos within x and is the total number of sub-halos inside -R200 • 



3.1. Point sources 

Under the assumption that the lenses do not overlap along the line of sight, the optical depth 
T represents the fraction of a given patch of the sky that is covered by regions in which a point 
source will be lensed. In the limit of small r, the optical depth can directly be used as an 
estimate of the lensing probability. 
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Figure 2. Optical depth for 
a point source at Zg = 1 as 
a function of projected radius. 
Here we assume a host halo 
at zi = 0.5 of mass Mhost = 
1.8 X 10^^ Mq and sub-halos 
in the mass range 4 x 10® - 
lO^'' Mq. The peak has been 
cut with respect to a maximum 
magnification factor of 50 at the 
Einstein radius of the host halo. 
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with the area of a patch on the sky S, the projected radius from the center of the host halo ^ 
and the external convergence k and shear 7 induced by the host halo. Here ciens denotes the 
cross section for a single sub-halo lens within the potential of the host halo. We use the analytic 
expression derived in [27] to compute the lensing cross section cjiens foi" SIS sub-halos. 

As can be seen in Fig. [21 the optical depth for a typical scenario does not exceed a value of 
4 X 10^'^ at any radius. This is contrary to prior claims where the optical depth for sub-halo 
lensing had been estimated to be several orders of magnitude higher [9]. This disagreement can 
be traced to an error in the sub-halo mass function adopted by the previous study. 

3.2. Extended sources 

According to |2SJ , it should be possible to detect lensing effects from dark matter substructures 
when observing extended sources resolved at scales smaller than the Einstein radii of the sub- 
halos. Furthermore, one should even be able to put constraints on the internal density profiles of 
the lensing sub-halos. This technique may already become observationally feasible with ALMA 
[15] or future space- VLBI missions like VSOP-2 [29] and thus constitute a major step forward in 
the study of dark halo substructures. Here, we estimate the probability of sub-halo lensing for 
e.g. a quasar in the radio regime by computing the average number of sub-halos that lie within 
the region of the host halo covering the source. 

In Fig. [3l the expected number of sub-halos covering an extended source at = 1 as a 
function of projected radius from the host halo lens center is shown for several source radii rg, 
ranging from 10 pc to lOkpc. This can be compared to the virial radius -R200 ~260 kpc for 
the host galaxy at zi = 0.5. It becomes clear that for sufficient source size (>lkpc) there is a 
good probability for the source image to be affected by sub-halo lensing, not only close to the 
Einstein radius of the host halo but even at a rather large projected distance from the host halo 
lens center. For a source with Vs = 1 kpc, one would expect at least one intervening sub-halo per 
10 observed systems with a maximum projected distance of 10" between the foreground galaxy 
and the source. For rg = 10 kpc, this number increases to approximately 10 sub-halos projected 
on the source out to a distance of 10" from the host galaxy. 
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Figure 3. Average number of 
substructures covering an extended 
source at = 1 as a function of 
projected radius for different source 
sizes. We assume a host at zi = 
0.5 with Mhost = 1-8 X lO^^M© and 
sub-halos in the range 4 x 10^ - 
10^^ Mq. We plot our results for 
a source with radius rg = 10 pc 
(dotted line), 100 pc (dash-dotted 
line), 1 kpc (dashed line) and lOkpc 
(solid line), respectively. The peak 
has been cut with respect to a 
maximum magnification factor of 30 
at the Einstein radius of the host. 
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Figure 4. Average number of 
substructures covering a source with 
radius = 1 kpc at a projected 
radius of R200/2 (ranging from 22"- 
7" for zi = 0.5 - 2) as a function of 
lens redshift, zi, and source redshift, 
Zg. We assume a host halo of mass 
M = 1.8 X IO^^Mq and sub-halos in 
the mass range 4 x 10^ - lO^^M©. 



We also explore the dependence of these results on the lens and source redshifts, zi and 
Zg respectively. The probability for a quasar to be aligned close to the Einstein radius of a 
foreground galaxy, where its magnification will boost the expected number of intervening sub- 
halo lenses, is low even when considering magnification bias. However, it is expected that there 
is a large number of galaxy-quasar pairs with a projected separation smaller than the virial 
radius of the host halo R200 (typically up to a few tens of arc seconds). In Fig. [H the average 
number of intervening sub-halos at i?2oo/2 for a source size of rg = 1 kpc is shown as a function 
of zi and Zg. The probability for lensing by sub-halos increases by up to one order of magnitude 
for high-z objects. 

Since high-mass host halos typically possess a larger fraction of their mass clustered in 
substructure, this will also boost the expected number of sub-halos projected on the source. 
Furthermore, high-mass host halos possess larger virial radii -R20O) allowing for larger projected 
distances between the host galaxy and the background source (see [30] for details). 



4. Discussion and Conclusions 

Our results indicate that the detection of sub-halos through gravitational image-splitting is 
likely to be considerably more challenging than suggested in previous studies, due to the 
smaller image separations predicted for sub-halo density profiles more realistic than the SIS 
models often adopted. In fact, no currently planned telescope will be able to resolve the image 
separations produced by sub-halos with density profiles of the type suggested by the most 
realistic simulations recently available (H03 & K04). If the sub-halos would have steeper central 
density slopes (e.g., of M99 type), these would give rise to image separations that could be 
resolved even with existing telescopes. A recent simulation with the highest resolution currently 
available [31], obtained central density profiles of the form p{r) oc r~'^ with central density slopes 
a ~ 1.24 (i.e., between the unfavourable NFW with a = 1 and the optimistic M99 profile with 
a = 1.5). If strong lensing effects produced by sub-halos which such profiles could be observed 
with any present or future telescope, is something we will be looking into in the near future. 

Despite the somewhat bleak detection prospects presented here, there are at least two effects 
that can potentially improve the detectability of image splitting by sub-halos: baryon cooling 
and the presence of intermediate mass black holes. 

Baryon cooling would cause the sub-halo (or its progenitor halo) to contract (e.g., [32t [33 tlMl 
I35j). thereby increasing the central density and boosting the image separation. It is, however, not 
clear how strong this effect is likely to be, since this depends on the details of the mechanism 
that prevents the dark sub-halo from forming stars. Baryon cooling is usually assumed to 
be associated with star formation, and most attempts to explain the lack of bright sub-halos 
therefore propose that the baryons have either been lost early in the history of the Universe, 
or are prevented from cooling by the ultraviolet background provided by re-ionization [36^ [37] . 
Neither mechanism is likely to result in any significant halo contraction. Nonetheless, claims 
of dark baryons in the form of cold gas in galactic disks have been made [38], implying that 
there may be some route for gas to cool without forming stars or being detected by the usual 
H2 tracers |39[ I40j. One may therefore speculate that there could be alternative scenarios, in 
which sub-halos are kept dark even though baryon cooling has taken place. 

Intermediate mass black holes (IMBHs, see e.g., |4 H \42 \ H3]) with masses of ~ 10^-10^ Mq 
would also boost the image separations, if present in the centres of dark sub-halos. This could 
for instance be the case if the empirical relations between the mass of a super-massive black hole 
and its dark matter halo [H] would extend into the dwarf-galaxy mass range. Even in the case of 
K = 7 = 0, an IMBH of mass 10"^ Mq, would give an image separation of « 4 x 10 ^" for z\ = 0.5 
and Zs = 2.0. This is sufficient to allow VSOP-2 to resolve the image splitting, regardless of the 
density profile of the sub-halo hosting the IMBH. Of course, even if the Ferrarese relations would 
hold for luminous dwarf galaxies, they do not necessarily do so for dark ones, since this depends 
on the formation details of IMBHs. Moreover, a population of halo IMBHs not associated 
with sub-halos could form an undesired background of milli-lensing events that would obfuscate 
attempts to study sub-halos through image-splitting effects. 

We have shown that the optical depth r for sub-halo lensing of point sources is lower than 
previously predicted. Even for the most favorable conditions it barely exceeds 0.01. We conclude 
that it is currently not feasible to use this technique to search for strong lensing signatures in 
point sources as e.g. quasars in the optical. 

If one instead targets extended sources, such as quasars in the radio wavelength regime, there 
is a high probability for sub-halo lensing of sources of sufficient size. For source sizes > 1 kpc, 
this is valid even at rather large projected distance of the source to the host halo center. This 
allows for a different search strategy than those previously proposed. Instead of only targeting 
multiply-imaged quasar systems, even quasar-galaxy pairs with a separation of several tens of 
arc seconds should show effects of strong lensing by substructures in the lens galaxy halo. 

However, these effects will strongly depend on the sub-halo mass and density profile. Even for 



Figure 5. Average number of 
substructures covering a source with 
Tg = 1 kpc at Zs = 1 as a 
function of projected radius for 
different minimum sub-halo masses 
Msub_min- For a host halo at z\ = 
0.5 with Mhost = 1-8 X lO^^^©, 
we vary Msub_min with 4 x 10^ Mq 
(dotted hue), IQ^Mq (dashed hue) 
and IQ^Mq (solid line), respectively. 
The peak has been cut with respect 
to a maximum magnification factor 
of 30 at the Einstein radius of the 
host halo. 

the most favorable density profiles, there will be minimum masses for which the image separation 
drops below the resolution of any present or planned observational facility. E.g., for the currently 
available EVN, the minimum sub-halo mass which could be resolved under optimal conditions 
is approximately 4 x 10^ Mq in case of an SIS profile but increases to 3 X lO^M© for an M99 
profile. Figure [S] shows how the expected number of intervening substructures N depends on 
the minimum sub-halo mass Msub_min- For a host halo of mass M = 1.8 x lO^^M© at z\ = 
0.5 and a source of size rg = Ikpc at = 1, we compute as a function of projected radius 
from the host halo center for minimum sub-halo mass -/Wsub_min = 4 x IO^Mq, 10^ and 
10^ Mq, respectively. Since the sub-halo mass function predicts about equal mass within each 
logarithmic mass interval, most sub-halos will be of low mass (n(Msub) oc Thus, the 

expected number of intervening substructures is very sensitive to the minimum sub-halo mass 
that can be resolved. 

Therefore, angular resolution will be crucial when attempting to detect CDM substructure 
via its lensing effects on background quasars and sub-milliarcsecond-resolution facilities will be 
required. We assess that the prospects for such a detection with the upcoming VSOP-2 satellite 
are good, with the condition that the ground array should be extensive enough. This is not a 
problem especially at 8.4 GHz where a number of large antennas are available from the geodetic 
VLBI network. However, one must be careful not to confuse internal structures found in quasars 
observed at radio wavelengths with sub-halo lensing signals. Quasars already macro-lensed on 
arc second scales can be used to test that it is possible to distinguish between the two, since 
internal structures should be mapped in all of the macro-images while sub-halo lensing will only 
affect one of the images. Taking the above into account, one should be able to use this technique 
to put constraints on dark matter sub-halos predicted by simulations in the near future. 
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